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Summary 


The  Air  For  ?e  has  a continuing  need  to  determine  the  response  of  electronics 
to  weapons,  space,  and  reactor  radiations.  Many  electronic  devices  consist  of  a 
combination  of  metallic  and  dielectric  elements.  The  response  of  metal-dielectric 
devices  to  radiation  is  not  well  understood.  Experimental  studies  indicate  that 
such  device  response  cannot  yet  be  predicted.  Problems  exist  with  transmission 
lines.  MOS-type  transistors,  capacitors,  large-scale  integrated  circuits,  and 
compton  diodes.  The  device  response  to  ’ adiation  has  sometimes  been  unpredict- 
ably  time  dependent.  Similar  devices  do  not  always  respond  in  similar  ways. 

A mechanism  which  may  help  to  explain  the  above  mentioned  difficulties  is 
developed  here.  Studies  on  high  energy  electron  transport  under  Air  Force  Proj- 
ect 5621  and  under  Defense  Nuclear  Agency  Subtask  Z99QAXTA040  and  by  the 
AEC  have  provided  data  which  allows  one  to  discuss  this  "new"  mechanism  quan- 
titatively. It  is  hoped  that  the  mechanism  developed  in  this  report  can  be  used  to 
help  clear  up  some  of  the  difficulties  in  predicting  device  response. 

High  energy  radiation  induces  electrical  currents  (sometimes  called  replace- 
ment currents)  in  metallic  elements  when  the  metallic  elements  are  adjacent  to 
dielectrics.  A cause  of  these  currents  is  known  to  be  the  radiation-induced  motion 
of  electrons  in  the  dielectrics.  The  high  energy  electron  transport  in  the  dielec- 
trics usually  produces  time  varying  electric  fields.  These  electric  fields  cause 
time-dependent  "conduction"  currents  in  the  dielectric.  Thus  the  currents  induced 
in  metallic  elements  are  time  dependent.  A mathematical  model  is  presented  to 
determine  the  time-dependent  currents  and  electric  fields.  Results  are  illustrated 
for  X-ray  radiations. 
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Radiation  Induced  Electrical  Current  and 
Voltage  in  Dielectric  Structures 


1.  INTRODUCTION 

The  purpose  of  this  work  is  to  theoretically  investigate  the  following  effect. 
High  eneigv  radiations  produce  electron  and  photo-compton-electron  currents  in 
all  materials.  Electron  currents  flowing  in  a dielectric  will  induce  electrical 
currents  on  adjacent  conductors.  Usually,  the  dielectric  currents  are  divergent 
and  will  deposit  electrical  charge  in  the  dielectric.;  This  charge  deposition  pro- 
duces electrical  fields  and  "conduction"  currents  in  the  dielectric,  and  thereby 
induces  further  currents  on  nearby  conductors.  We  develop  a calculational  model 
to  predict  the  induced  currents  to  the  conductors.  The  work  is  limited  to  a simple 
"one  dimensonal"  geometry,  which  can  be  compared  to  controlled  experiments 
and  which  can  be  understood  without  advanced  mathematical  techniques. 

The  simplest  arrangement  of  dielectrics  and  conductors  which  can  be  analyzed 
and  experimentally  controlled,  in  one  dimension,  is  shown  in  Figure  1.  All  work 
reported  here  is  limited  to  this  geometry.  Extension  to  more  complicated  geo- 
metry or  higher  order  dimens iorality  results  m greatly  increased  mathematical 
and  experimental  complexity.  Itemoval  of  any  of  the  elements  in  Figure  1 results 
in  an  arrangement  winch  cannot  be  useful  in  a real  device  application. 

(Received  for  publication  21  November  1974) 
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E3  CONDUCTOR 
□ DIELECTRIC 

(T)  ,-AlPERL  ;.1*TER,  Arrow  shorn  positive  current  direction, 
ant  surtece  -*v,rg«  density,  coulombs/meter^ 
pat  dielectrn  rherge  density,  coulombs/meter^ 

Radiation  travels  it.  positive  X direction. 

O.a.b.c  are  the  X posi'.ons  of  the  dielectric  - metal  interfaces 
E is  electric  field,  volts/meter. 


Figure  1.  The  Simpltrt  Metal  Dielectric 
Structure 


2.  ELECTROSTATICS 

Consider  the  conditions  prior  to  t'me  t = 0 sec.  Vj  and  V2  have  been  applied 
and  quiescent  «r0<t  = 0),  0^=0),  Dgit-- 0),  and  <73(t=G)  have  been  established.  Mo 
current  flows  in  the  meter,  I , since  the  dielectric  polarization  has  equilibrated 
and  since  the  dielectrics  are  not  "leaky".  At  t !=  0 sec,  we  apply  the  radiation 
which  passes  through  the  material,  creating  an  electron  current  density  J(x,  tS  0) 
in  the  x direction. 

This  current  density  can  be  described  by  sheets  of  charge,  p(x,  t)dx  moving 
with  some  velocity  along  the  x axis.  Every  charge  sheet  is  the  source  of  induced 
electric  fields  and  Ep.  The  line  integral  sum  equation 

El  • dl  + p Er  • dl  = 0,  (11 

0 x 

is  always  satisfied.  A surface  charge  density,  Ad,  is  induced  at  the  conductor 
surfaces  in  response  to  the  sheet  of  charge  such  that  the  normal  component  of  the 
electric  field  vanishes  in  the  conductor,  and  we  have 


10 


+ Affft(t)  _ - Aff.(t) 

EL(t)  = x ai.d  ER(t)  = x . (2) 

where  x is  the  unit  vector  and  £ is  the  permittivity  of  the  dielectric  in  F/m.  By 
applying  a gaussian  surface  arocrd  the  entire  system  we  can  show  that 

Ao0(t)  + - -p(x.t)dx.  (3) 

Combining  Eqs.  (1),  (2),  and  (3)  results  in,  for  0 s x < a, 

Ao0(t)  = p(x,t)dx  and  t'  = -|[p{x,t)dx  (4) 

and  by  a similar  process  we  get,  for  b * x < c , 

Aa2(t)  = ^-|p<x,t)dx  and  ikCgjt)  = ^|p(x,t)dx.  (5) 


As  p(x)  is  distributed  over  all  x , «e  can  , -**  p over  all  x to  obtain  :he  total  in 
duced  conductor  surface  charge.  For  example. 


a 

EAffjit)  = y - — p(x,  t)dx 
0 

Ojit)  * o.(t=0)  + \ - ^p(x.t)dx  , 

0 


(6) 


and  similarly, 

c 

EA°2(t)  = J fffp(*.t)dx 

0 

c 

»2(t)  a ff2^=0^  + J fr§  p<*.  t)dx  • 

b 


3.  ELECTRODYNAMICS  4 

The  experimentally  measurable  parameter,  I , is  composed  of  the  currents 
required  to  establish  Oj(t)  and  and  the  currents  actually  injected  across 

the  metal-dielectric  interface.  Thus,  in  our  one  dimensional  case. 
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m>  = — ar~  a r~ 

where  J is  in  units  of  A/m  . 
From  Eq.  (6)  we  get 


do  ^t) 


d02it) 


+ J(a.t)  - J(b,t) 


do,  (t) 


IT  if 

It  ~ " a J 


ap(x,t> 

x — — — dx  . 


The  continuity  equation,  - = V • J , in  the  one  dimensional  geometry  of  this 

problem  gives 


do.  (t)  . oa  dJ'x.t) 

TT = — \ X — ^ dX. 

it  a J dx 

0 


Integration  by  parts  results  in 


do,  (t) 


a 

= J(a,t)  - ^ ^ J(x,t)dx  , 


and  simi..  ly, 


do2(t) 


Jib,  t)  + y Jix,t)dx  . 


Substituting  into  Eq.(8)  above,  we  get 


- i ^ J(x,t)dx  - ^ J(x,t)dx. 


Thus  we  have  shown  that  the  meter  current,  lit),  depends  entirely  cn  the  volume 
integral  of  cui  rent,  J(x,t)  in  the  dielectric.  In  cases  where  we  can  estab’ish  the 
current  in  the  dielectric,  under  radiation,  we  can  predict  the  meter  response. 
Calculations  and  measurements  of  radiation -induced  currents  are  being  done  else- 
where and  will  be  referenced  later  in  thia  report.  We  can  assume  for  the  moment, 
that  J(x,t-0)  is  known. 
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A new  problem  arises  since  we  need  to  know  J(x, t > 0).  The  application  of 
radiation  to  solids  produces  conduction  electrons  (as  well  as  "holes")  which  can 
turn  a dielectric  into  a semiconductor.  Usually,  J(x,t  = 0)  is  such  that  V • J 4 0, 
and  charge  will  be  deposited  at  a rate  given  by  = - V • J(x,  t).  This  charge 

deposition  creates  a changing  electric  field  in  the  dielectric,  and  the  radiation- 
inducea  conduction  electrons  move  in  response  to  the  new  electric  field  strength. 
Thus,  for  t > 0 we  need  to  include  this  conduction  current  in  th?  determination  of 
J(x,t  > 0). 

We  shall  calculate  this  conduction  current  using  the  following  steps: 

(a'  Let  n = 0,  where  n is  a running  integer. 

(b)  Calculate  the  electric  field  produced  by  the  charge  deposited  through  a 
"short"  time.  At. 

(c)  Determine  the  conduction  current,  jc<x,  t),  produced  by  the  electric  field 
at  the  end  of  the  time  interval,  t = (n+l)At. 

(d)  Assume  that  at  t = (n+l'M,  a new  curr*>"t.  is  flowing  defined  by 

J(x,  t)  = J(x)  + J (x,  t).  We  will  discuss  J (x,t)  and  J(x)  later, 
c c 

(e)  Step  the  integer  so  that  n+  1 = n and  return  to  (b)  above. 

In  this  way  we  are  able  to  determine  the  time  dependence  of  J. 

We  determine  the  steps  in  the  following  ways. 

To  calculate  the  electric  field,  (b),  refer  to  Figure  2.  We  calculate  E(x) 
using  Gauss'  law, 

£ E • n dS  = j*  jElll  dv  , 
v 

applied  to  the  "gaussian  surface".  Since  the  apparatus  extends  infinitely  in  the  y 
and  z directions,  symmetry  considerations  imply  that  E • n must  be  zero  every- 
where on  the  upper  and  lower  sides  of  the  "gaussian  surface".  E • n is  also  zero 
everywhere  in  the  conductors. 

Thus  we  have 


f 


E • n dS  = AE(x) 


where  A is  the  cross  section  area  of  the  "gaussian  surface"  at  x. 

At  depth  Xj  in  the  dielectric,  E is  produced  by  charges  enclosed  in  the  two 
gaussian  surfaces  shewn  in  Figure  2.  The  total  field  at  x;  is  then  given  by: 


E(x  ,t)~T 

x t 


Xj  a 

J p(x,  t)dx  - ^ p(x. t)dx  + ffg(t)  - ffj(t) 


+ Eext(xl> 


(14) 
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Figure  2.  Calculation  of  Electric  Field  F rom  Deposited 
Charge  Distribution 


where  Eext<*j)  is  the  electric  field  provided  by  external  voltages  as  shown  in 
Figure  1.  Note  that  Ee^(x)  a E(x.  0). 

The  function  p(x.t)  can  be  approximately  determined  from  the  equation  of 
continuity  by 

p(x.t)  = -V- J(x.t')dt'  rp(x.t-At)  - At  [v • J(x, t - At)J . (15) 

t-At 

In  this  formalism  we  require  that  J(x,t')  not  vary  significantly  (less  than  1 percent 
or  so)  during  the  time  span  t - At  < t'  < t.  We  determine  [see  Eq.  (11)]  after  a 
short  time  interval.  At,  by: 


<r,(t)  = aj(t-At)  + | jj(a, t')  ~ i y J(x,t')dxjdt' 

Oj(t)  = o^t-At)  + At  £j(a,t  - At)  - £ J J(x,t  - At)dxj, 

and  similar  equations  give  us  «tq,  or2>  and  ffg  at  time  t.  Thus,  for  example, 
the  field  E(Xj)attime  At  is  given  by  Eqs.  (14-17)  as 

rx  1 

C E(Xj , At)  = € E(Xj  , 0)  + \ p(x,  0)dx  - At  ^J<Xj , 0)  - J(0 , 0)J 
ra 

- \ p(x,  0)dx  + At[j(a,  0)  - J(x1(0)]  4 oQ(At)  -ffj(At) 


(16) 


(17) 


(18) 
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where  E(Xj , 0)  s e ^ evaluated  at  Xj . We  will  generally  assume  that  prior  to 
t - 0,  the  dielectrics  contain  no  charge;  p(x,  0)  = 0.  It  is  only  necessary  to  know 
E(x,  0)  and  J(x,  O'  at  time  t = 0 to  determine  E(x,  At)  as  long  as  At  is  sufficiently 
short. 

Now  we  must  determine  (c),  .Ic(x,t).  In  general,  for  constant  intensity  radiation, 

J(x,  t)  s J(x)+Jc(x,t)  (19) 

where  J is  a conduction  current  which  is  driven  by  the  electric  field  and  J(x)  is 
c 

the  field  independent  radiation-induced  high  energy  electron  current.  We  define 
Jc(x,t)  h[«0+  «RD(x)]E(x.t)  (20) 

where  D(x)  is  the  radiation  dose  rate  in  rad /sec,  ZJq  is  the  normal  conductivity 
of  the  dielectric,  and  ISR  is  the  coefficient  of  radiation-induced  conductivity. 

D(x),  like  J(x),  is  determined  by  methods  in  the  literature  which  will  be  refer- 
enced later.  Published  experiments  indicate  that  radiation-induced  conductivity 
is  proportional  to  dose  rate  in  many  materials. 

We  are  able  to  determine  I(t  = At)  from  the  above  equations  by  evaluating 

J(x,  At)  = ,I(x)  + E(x,  At)[«0  + UrD(x)]  ,21) 


where  E(x,  At)  is  obtained  from  Eq.  (14)  above  for  0 < x < a and  from  a similar 
equation  for  b < x < c , and  plugging  J(x,  At)  into  Eq.  (13)  to  get 


I(At) 


i f J(x;  At)dx  - —t-  f J(x,  At)dx. 

a ■'n  c - D 


(22) 


The  cycle  (d,  e)  can  now  be  completed.  We  have  calculated  I(t=  0)  and  I(t=  At). 
To  calculate  I(t=  2 At)  we  have  only  to  go  through  the  calculation  again,  being 
careful  to  replace  t = 0 by  t = At  and  also  t = At  by  t = 2 At;  and  by  iterating  we 
can  obtain  I(nAt)  for  any  reasonable  number  n. 

A computer  program  has  been  written  to  do  this  calculation  by  the  iterative 
technique  given  above.  Thus  we  can  determine  I(t),  J(x,t),  E(y,t)  for  any  type 
of  radiation  which  produces  a known  J(x,  0)  and  D(x,  0)  in  this  one  dimensional 
geometry. 

A flow  chart  for  the  computer  program  is  given  in  Figure  3. 
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USE  EQS.I18) 
TO  DETERMINE 
E(x.  0 


USE  EQ.  (21 ) 
TO  DETERMINE 


USE  ECM13) 
TO  DETERMINE 
I (tl 


t = NAt 


Eq.  (13)  I («)  * j-  f J(x,  tlrix  - -1-  fj(x.  tldx 

(a0M‘  CT0(At(N-1l)+  At  [,,0  At(N-ll)  4 J-  j J |x.  (N-1)  At  )dxj 

^ ^ ((N-DA«)+  At  [j(..  At  (N-ll)  - I f J | x.  At  (N-1)  j dxl 

* 'o  J 

CT2(t)=  a2  ( (N-1)  At)+  At  |^-J{  b,At  (N-l) ) + ± f J (x.  At  (N-1))  dxl 

a3  (*»  * a3  ( <N-D  At)+  At  [j  (c.At(N-D)  --L  J J { x.  At  (N-1)  )dxj 

a0  ,01  * °i  <0*  1 °2  ,0)  * °3  ,0)  * 0 


Eq  (18) 


i x < i , E (x,  t)  - 1 |(7q  (t)  - •Tj  (t)  + A t [j  (i,  m A t)  - 2J  (x,  m A t)  4 J (0,  m A till  + E (x  ( 
L m 0 J 

< x i e;E(x,t)»i  <t)  _ »3  (,)  + At  £ [j  (c,  m A t)  - 2J  (x.  m A t)  + J (b,  m A t)]j+  E (x.  ( 


Eq.  (21)  Jlx.tl  ■=  J (x)  4 E(x.t)ju0  4 UR  D (xlj 


Figure  3.  Computer  Program  Flow  Chart 
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The  analysis  presented  is  a straightforward  application  of  electrostatic  equa- 
tions to  a one-dimensional  linear  dielectric  material.  Similar  analyses1  have 
been  published  for  various  applications  (such  as  vacuum  tubes,  semi-conductors 
and  others)  and  can  be  found  in  many  textbooks  on  applied  electrodynamics.  We 
have  included  a derivation  of  the  equations  to  make  clear  the  exact  conditions  and 
assumptions  with  which  we  are  working. 

The  equations  are  valid  under  the  following  conditions: 

(a)  £ is  constant  throughout  the  dielectric. 

(b)  The  dielectric  is  in  reasonable  contact  with  the  metal  plates  during 
irradiation.  Radiation-induced,  low-energy  electron  emission  effectively  short 
circuits  any  gaps  between  materials. 

(c)  The  electron  flux  directly  generated  by  the  high  energy  radiation  is  not 
significantly  perturbed  by  the  electrostatic  fielda. 

(d)  J and  E are  functions  of  only  x and  t. 

(e)  The  radiation- induced  conductivity  must  be  known  and  must  vary  slowly 
with  E;  that  is,  no  avalanching  or  breaking  down  effects  are  considered. 

Electron  currents  generated  within  the  dielectric  by  the  radiation  create  the 
new  ana  peculiar  twist  shown  by  these  equations.  Moving  electrons  can  originate 
and  "die"  entirely  within  the  dielectric.  Generally,  the  energy  distribution  of 
these  electrons  is  such  that  the  electric  fields  do  not  substantially  alter  the  tra- 
jectory oi  the  vast  majority  of  the  high  energy  electrons.  The  radiation-induced, 
high-energy  currents  .I(x)  can  be  calculated  from  the  photoelectric,  compton,  and 
high  energy  electron  transport  effects.  Radiation  simultaneously  excites  low 
energy  electrons  which  are  strongly  affected  by  the  electrostatic  fields. 


4.  DETERMINATION  OF  THE  PHOTO-COMPTON  CURRENT,  J(x) 

4.1  Hiitorical  Background 

A historical  perspective  aids  one  in  understanding  the  approach  in  this  report. 
High  energy  radiation  transport  has  been  a large  important  field,  both  economically 
and  politically,  for  30  years.  The  major  radiation  interactions  can  be  quantitatively 
discussed  such  that  many  engineering  solutions  can  be  deduced  from  physical 
laws.  Many  important  economic  and  military  problems  have  been  heavily  studied. 
Specifically,  the  design  of  fission  reactions  to  optimize  output,  the  design  of 

~ M;  rx,  K.  D.  (1973)  IEEE  Trans.  Nuc.  Sci.  NS20(No.  6h64-71.  (In  this  paper, 
equations  are  derived  for  a t rans  mission^  ine  under  conditions  which  could 
be  caused  by  radiation.)  Also  se?  Birdsall,  Charles  K. , and  Bridges, 
William  B.  (1966)  Electron  Dynamics  of  Diode  Regions,  Academic  Press 
Inc. . New  York.  N."Y^ 
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shields  to  atteuuulc  penetrating  radiations  from  reactions,  the  calculation  of  apace 
and  atmospheric  transmission  of  nuclear  radiations,  the  Interaction  of  radiations 
with  biological  tissue,  the  use  of  radiation  processing  to  create  better  materials, 
and  the  displacement  of  atoms  in  solids  are  alt  problems  of  strong  continuing 
Interest. 

A great  deal  of  effort  hue  produced  useful  quantitative  Information  which  can 
be  quickly  applied  to  many  specific  engineering  problems.  Most  of  this  informa- 
tion Involves  an  energy  flux  or  an  energy  transport  since  many  of  the  above  men- 
tioned effects  can  be  approximately  related  to  the  rndiatlon  energy  density. 

Nearly  always,  the  energy  transport  is  nccomponled  by  a current  flow  of  high 
energy  electrons.  It  Is  rare  that  the  current  flow,  1 , is  responsible  for  tissue 
damage,  reactor  shielding  problems,  or  ntomic  displacements!  the  energy  den- 
sity term  correlates  more  strongly  with  these  effects.  But  we  now  have  computa- 
tional techniques  for  determining  the  energy  density  for  many  kinds  of  indlotion 
In  mnny  arrangements  of  matter. 

From  the  computational  techniques  developed  for  determining  the  radiation 
energy  density,  we  can  determine  the  current  flow  .l(x>  in  some  geometries.  It 
is  not  a simple  problem.  The  present  state  of  the  nrt  Is  best  displayed  In  Refer- 
ence  2 and  In  the  references  contained  in  Reference  2.  In  this  reference  a cur- 
rent .Hx'  is  calculated  for  monoenergetlc  monodlrectlonnl  photons  in  homogeneous 
material.  .I(x)  Is  tabulated  for  25  elements  and  10  Insulating  compounds  under 
photon  radiation  of  energies  from  0.01  MeV  to  20.0  MeV.  However,  the  restric- 
tion to  homogeneous  materials  is  too  severe  for  most  applications.  An  example 
of  application  to  a nearly  homogeneous  apparatus  Is  gB'cm  In  Reference  3. 

The  situation  for  non-homogeneous  material  has  beer,  generally  "Ignored"  by 
analytical  workers  since  the  problem  appears  to  be  very  complex.  Occasionally, 

apparently  motivated  by  some  application,  experimentalists  hove  made  measure- 

4 5 6 

ments  of  .l(x>  at  on  Interface  between  two  metnls.  ' 

A large  number  of  measurements  under  non-rigorous  conditions  were  made 
of  the  current  flowing  to  metal  plates  Imbedded  in  a dielectric,  thus  attempting 
to  use  the  device  as  * radiation  detector;  the  exact  dependence  of  7(x,y,  z,t)  wns 
not  known.  For  information  on  these  metal-dlelectrlc  experiments,  see  Refer- 
ence 7 and  especially  see  the  references  listed  extensively  therein. 

2.  MacCallum,  Crawford  J,  and  Dellin,  Theodore  A.  (1973)  .1.  Appl.  Phys.  44 

(No. 41:1878-84.  * 

3.  Brumley,  F.  B.,  et  al  (1073)  IEEE  Trans,  Nuc,  Scl,  NS-20(No.  81:48-57. 

4.  Hess,  Bernhard  (1050)  Zelt.  Ang,  Physlk,  Band  11  Heft  j^:449-53, 

5.  Rnab,  Bernard  (1083)  Nucleonics  21(No.  21:46-7. 

8.  Ebert,  P..F,  and  Lauzon,  A.  F,  (1087)  Rev.  Scl,  Inst.  3B(No.  121:1747-52. 

7.  Fewell,  T.  R.  (1072)  Compton  Diodes;  Theory  and  Development  for  Radiation 
Detectors,  Sandia  Labs,  Albuquerque,  N.  M. , Report  S>C-i)R-7zoi  18. 
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Approximate  analytical  techniques  have  been  devised  to  make  rough  calcula- 
tions of  the  currents  near  an  interface.  Many  works  have  assumed  that  an  elec- 
tron is  initially  generated  by  the  photoelectric,  compton,  and  Auger  effects  and 
then  travels  in  a straight  line  until  it  stops  one  range  away.  Such  a calculation 
can  only  give  order-of-magnitude  estimates  for  the  spatial  average  current  and 
can  often  give  the  wrong  sign  of  I.  Near  an  interface,  the  multiple  scattering  of 
electrons  and  the  discontinuous  change  in  the  high  energy  electron  generating  func- 
tion can  actually  reverse  the  direction  of  J(x)  over  a finite  spatial  interval.  The 
equations  developed  earlier  indicate  that  we  need  to  know  div  J(x,  t)  to  determine 
the  time-dependent  currents  and  fields  in  a dielectric;  this  requires  relatively 
accurate  knowledge  of  the  spatial  dependence  of  J(x). 

O 

Work  is  progressing  at  this  laboratory  and  elsewhere  to  measure  the  spatial 
dependence  of  .1  in  planar  goemetry  for  non- homogeneous  slabs.  Calculational 
schemes  are  being  developed  to  determine  J(x,  y,  z.t)  for  more  complex  geometry. 
The  planar  experiments  will  be  used  to  test  or  validate  the  computations.  When 
the  computations  properly  predict  the  planar  experiments,  they  may  be  used  to 
determine  the  time-dependent  current  for  radiation  spectra  and  more  complex 
geometries  which  are  of  interest  to  the  Air  Force. 

4.2  Approximate  Function*  for  J(x) 

Short  of  having  exact  functions  for  J(x)  determined  by  experiment  and  compu- 
tation. we  give  approximate  J (x)  for  a few  geometrical  cases.  Similar  approxima- 
tions can  be  found  described  in  Reference  9.  In  these  curves  we  neglect  the  rela- 
tively small  effect  due  to  photon  attenuation. 

We  consider  1.25  MeV  gamma  rays  and  0.2  MeV  gamma  rays  incident  m the 
normal  direction.  The  shapes  presented  need  not  be  accurate  since  our  intent  is 
not  to  make  exact  predictions.  The  curve  shapes  are  only  indicative  of  the  kinds 
of  effects  seen  with  gamma  rays  at  these  energies.  The  curves  are  described  in 
Figure  4. 

In  Figure  4,  the  number  z is  numerically  equal  to  the  depth  in  the  dielectric 
expressed  in  meters.  It  is  estimated  that  J(x)  in  these  functions  is  everywhere 
within  50  percent  of  being  correct  for  the  given  y-ray  in  dielectric  materials 
with  effective  Z 4 18  near  metals  with  effective  Z £ 40.  Most  of  the  fits  are 
within  ±20  percent  of  the  true  data  over  the  majority  of  the  depth. 

The  data  in  Figure  4 applies  to  a single  metal-dielectric  interface.  We  can 
obtain  J(x)  for  the  apparatus  of  Figure  1 by  a superposition  of  forward  and 

8.  Kooi,  C.F.  and  Kusnezov..  N.  (1973)  IEEE  Trans.  Nuc.Sci.  NS-20(No.  6):97-104. 

9.  Chadsey,  William  L.  (1973)  Mcnte  Carlo  Analysis  of  X-Ray  and  y-Ray  Tran- 

sition Zone  Dose  and  Photo -T  inipton  Current,  AFcRl-tR-73-05?2. 
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Figure  4.  Analytical  Fits  to  Photo-Compton  Current  Data 

backward  terms  only  for  the  case  ir  which  the  dielectric  and  metal  are  each  at 
least  1/4  electron  range  thick.  Thinner  slabs  have  not  yet  been  analyzed  to  ac 
curately  determine  the  form  of  J(x).  We  will  therefore  always  work  with  slabs 
of  approximately  t electron  range  thickness. 
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J(x)  for  thinner  slabs  is  currently  being  experimentally  det  ir  mined  by  this 
author  and  analytical  v by  others.  The  thin  slab  arrangement  (less  than  one  range 
thick)  will  be  studied  i\'  the  near  future. 

J(x)  for  other  high  energy  radiations,  such  as  electrons,  is  available  but  is 
not  considered  in  this  report. 


5.  RADIATION-INDUCED  CONDUCTIVITY 


To  determine  the  conduction  currents  in  the  dielectric,  it  is  necessary  to 
know  the  radiation- induced  conductivity  everywhere  in  the  dielectric.  The  high 
energy  radiations  produce  a quasi-contmuous  energy  distribution  of  excited  elec- 
trons. The  distribution  extends  from  the  valence  band  to  the  quantized  energy  of 
the  incoming  radiation.  All  of  these  excited  electrons  respond  to  the  electric 
field  and  produce  a net  conductivity  current.  In  addition,  the  radiation  produces 
"holes"  which  can  similarly  respond  to  electric  fields.  The  problem  of  determin- 
ing the  radiation- induced  conductivity  from  the  excited  electron  energy  state  dis- 
tribution is  beyond  the  scope  of  this  report. 

We  shall  take  a simple  empirical  approach.  Induced  conductivity  has  been 
measured  in  thin  insulating  films  under  external  applied  bias.*®  Ti.e  photoconduc- 
tivity can  be  approximated  by  a linear  function 

Photoconductivity  = Up  D(x) 

where  Up  is  a constant  of  the  material  in  units  sec  ohm  1 meter  1 rad  *,  and 
D(x)  is  in  units  rad/sec. 

The  data  available  for  organic  dielectrics  (with  dielectric  constant  between  1 
and  10)  indicates  that 
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2.6  X 10  “ -r  £ U„  s 2.6  X 10"”  ~r -as= -? 

ohm  m rad  R ohni  m rad 


-14 


sec 


The  value  of  Up  for  any  one  material  has  varied  typically  a factor  of  10 
from  sample  to  sample  and  from  experiment  to  experiment.  The  causes  of  the 
variance  are  not  known.  Certainly,  low  energy  electron  transport  plays  an  im- 
portant role.  Experimental  errors  are  also  likely  to  have  caused  some  of  the 
discrepancies  in  the  value  of  Up  . 

The  work  in  this  report  does  not  depend  on  knowing  Up  accurately.  We 
shall  look  at  effects  throughout  the  probable  range  of  Up  given  above.  Up  is 

10.  Weingart,  I!'  C. , et  al  (1972)  IEEE  Trans.  Nuc.  Sci.  NSl9(No.  6):  15-22. 
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an  interesting  factor,  but  it  is  not  critical  to  an  understanding  of  the  qualitative 
features  of  this  work.  It  may  turn  out  that  the  experiments  proposed  here  could 
be  used  to  more  accurately  determine  0^  for  any  dielectric,  but  that  is  beyond 
the  sccpe  of  this  report. 

it  is  essential  that  ISp  be  effectively  const ar.<  throughout  any  one  experimental 
sample.  Without  this  homogeneity  all  the  foregoing  analysis  of  induced  conductivity 
is  inappropriate. 


6.  DEPTH-DOSE  DEPENDENCE,  D(x) 

The  energy  deposited  in  matter  by  radiation  has  been  extensively  studied. 
Recent  calculations  and  experiments  have  determined  dose  vs  depth  for  mono- 
direcaonal  radiation  normal  to  a planar  interface  of  two  materials.11  The  dose 
rate  D(x)  given  in  this  report  is  obtained  directly  or  by  analogy  from  published 
and  unpublished  works  of  other  people.  References  are  given  where  appropriate. 

By  varying  the  spectra  of  incident  radiation,  we  can  vary  the  dose  profile  sub- 
stantially. The  dose  profile  may  have  a strong  effect  on  the  current  flow  in  tbe 
apparatus  of  Figure  1.  We  shall  hypothesize  dose  profiles  which  may  or  may  not 
be  possible  in  order  to  determine  whether  any  unexpected  effects  are  likely.  It 
has  already  been  shown  that  the  dose  can  vary  a factor  of  30  or  more  with  low 
energy  X-rays  at  a high  Z - low  Z interface  over  a distance  short  relative  to  one 
radiation  attenuation  length.  In  princip’e,  b>  mixing  low  energy  (*  100  keV)  X-rays 
with  high  energy  (£  1 MeV)  X-rays,  it  is  possible  to  have  strongly  varying  dose 
profiles  and  weakly  varying  photo- oompton  current  J(x)  profiles. 


7.  RESULTS 

We  have  determined  I(t),  Eix.t),  and  J(x,t)  for  several  material  arrange- 
ments. The  radiation-induced  conductivity  coefficient,  , has  been  assumed  to 

- 16  ** 
be  10  sec/ohm  m rad,  since  this  appears  to  be  near  the  center  of  the  range  of 


11.  The  references  to  this  work  are  rather  diffuse.  One  current  starting  point 
for  is  'r.-rnation  is: 

Chadsey,  V/inicm  L.  <1973)  Monte  Carlo  Analysis  of  X-Ray  and  y-Ray  Tran- 
sition  Zone  Dose  and  Photo-Compton  Current,  AFCRL-Tll-', "3-0572,  espe- 
Hal'ly"  Refs.  1-4,  7.  S.  19. 

Another  nonoverlapping  source  is- 

Eisen,  H. , Rosenstein,  M. , and  Silverman.  J.  (1972)  Radiation  Research  52: 
429. 

And  another  is: 

Dutreix,  J.  and  Bernard,  M.  <1968)  BSophysik  4:302. 

The  radiation  shielding  community  has  extensive  dose  depth  data  for  homo- 
geneous materials,  and  the  radiation  therapy  field  has  studied  skin-air, 
and  tissue-bone  interface  profiles. 


22 


probable  values  for  organic  dielectrics.  The  permittivity  is  assumed  to  be 
€ = 2 X 10  11  F/m.  The  external  voltage  is  zero  in  all  cases. 

7. 1 Entirely  Homogeneous  Material 

Let  the  metal  and  dielectric  materials  be  of  the  same  effective  atomic  num- 
ber. The  photo-compton  current  profile,  J(x),  can  be  obtained  as  follows: 

(1)  Determine  the  photon  number  angle  en»rgy  spectra  vs  depth  in  the 
material  from  the  literature  on  X-ray  and  y-ray  attenuation  and  shielding. 

(2'  Determine  the  resulting  depth  dependent  photo-compton  current  aing, 
say.  Ref.  2 . 

The  dose-depth  dependence  can  also  be  obtained  from  the  extensive  literature  on 
radiation  shielding.  For  the  homogeneous  case  of  dimensions  given  below,  a 
simple  exponential  dependence, 

.Mx)  = .1(0)  e'^* 

D(x)  = D(0>e'4x 

is  approximately  correct.  The  approximation  is  sufficient  for  the  purpose  of  this 
report.  The  attenuation  coefficient,  p , is  given  in  meters  1 and  varies  strongly 
with  photon  energy. 

The  results  are  given  in  Figures  5 and  6. 

14  2 

Figure  5:  1.25  MeV  photons,  10  photons/m  -sec  on  ethylene. 

D(0)  5 rad /sec 

.1(0)  =-  -1.3  X 10'7  A/m2 
M = 6 m 1 . 

4 5 

E(x,t)  and  I(x,  t)  are  graphed  in  Figure  5 at  t = 0,  t = 10  , and  t^  10  sec. 

The  net  current  to  the  meter,  I(t)  = -7.4576  X 10  9 at  t 10® , and  I(t  = 0)  - 
-9  2 

-7.4582  x 10  A/m  . There  is  very  little  time  dependence  in  the  meter 
current. 

— 14  2 

Figure  6:  0.2  MeV  photons,  10  photons/m  -sec  on  ethylene. 

D(0)  = 1 rad/sec 
J(0)  = -4.5  X 10'9  A/m2 
(i  = 14  m * . 

The  net  current  to  the  meter,  I(t),  vane*  little  with  time;.  I(t  = 0)  = -.  4597  87  X 
10'10  A/m2  and  !<t  > 106)  = -.499786  X 10_1P  A/m2. 
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Figure  6.  0.2  MeV  Photons  on  Homogeneous 
Ethylene,  10^“*  Photons /m^ -sec 
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Since  the  material  is  homogeneous,  there  are  no  steep  curves  for  D(x)  or 
J(x)  and  the  data  is  relatively  simple.  These  graphs  are  included  so  that  com- 
parison of  magnitudes  may  be  made  between  the  homogeneous  and  non-homogene- 
ous  case.  Homogeneous  materials  of  effective  atomic  number  different  than  ethyl- 
ene give  qualitatively  similar  results. 

We  next  consider  situations  where  D(x)  and  J(x)  are  strongly  varying  functions. 

7.2  Center  Electrode  of  Differing  Atomic  Number 

Consider  the  apparatus  of  Figure  7.  We  can  vary  Zj  and  Zg  to  obtain  differ- 
ent effects.  As  always,  the  grounded  elements  must  be  good  conductors  and  the 
remaining  two  elements  are  dielectrics. 


z,  z,  z2  z,  z, 


Figure  7.  Apparatus  With  Center  Electrode  of  Differing 
Atomic  Number,  Z 

<1>  z2  » Zj 

14  2 

(a)  10  photons/m  -sec,  1.25MeV,  at  x = 0,  moving  in  positive  x 
direction.  Zt  — 5 and  Z„  = 29  sc,  that  the  materials  are  similar  to  ethylene 
and  copper.  Let  Z ^ have  a density,  p = 1.0  g/cm  and  a dielectric  constant 
€ = 2 X 10"11  F/m. 

The  functional  fits  for  J(x)  and  D(x)  are: 


0 * X £ a; 

J(x)  = 

[-  9.4  + 

1.6/ 

^(o.25  + 0.55  (0.005  -x)2  106) 

j 1.4  * 

© 

1 

00 

b ' x * c ; 

J(x)  = 

£-  9.4  + 

US/ 

^0.55  + (l.lx  -.Oil)2  106)] 

1.3  * 

CO 

< 

o 

m 


m 


25 


C £■  X * a ; 


D(x)  = 5 [l  + 1 /(l  + (0.005  -x)2  --  4 * 105)]  ||| 

b«x<c;  D(x)  = 4.8  [l  - 0.3  /(l  + (x  -0.01)2  * 3 * 105)]  |||  . 

The  'urve  fits  to  dose  depth  are  taken  from  the  data  of  Ref.  12  with  modification. 

The  dose  data  of  Ref.  12  near  depth  x = b (just  downstream  from  the  high  Z mate- 
rial) contains  a region  of  negative  slope.  It  is  assumed  that  this  region  of  Ref.  12 
data  is  due  to  low  energy  photons  (^400  keV)  causing  electron  photoemission 
from  the  high  Z into  the  low  Z material.  This  low  energy  photoemission  com- 
ponent of  the  dose  has  been  subtracted  out  to  obtain  dose  due  to  1.25  MeV  photons 

alone.  Monte  Carlo  calculations  of  dose  near  an  interface  strongly  support  this 

13 

subtraction  of  the  alleged  low  energy  component.  The  dose  depth  curves  of 
Ref.  13  are  also  similar  to  the  fits  given  above. 

The  results  for  Zg  » Zj  are  shown  in  Figure  8.  Two  especially  interesting 
results  obtain: 

g 

11  the  electric  fields  approach  10  V/m  which  could  cause 

dielectric  breakdown  in  some  materials,  and 

2)  the  meter  current  begins  positive  and  then  becomes  negative 

as  the  exposure  progresses  in  time. 

14  2 

(b)  200  keV  photons,  10  photons/m  -sec  incident  in  positive  direction 

13 

at  x = 0.  The  approximate  fits  to  the  available  data  for  J(x)  and  D(x)  can  be 
seen  in  Figure  9.  Note  that  J(x>  = J(x,  0). 

The  resulting  I(t),  Kx.t),  E(x,t)  are  plotted  in  Figure  9.  The  most  interest- 
ing features  are: 

g 

1)  The  maximum  electric  fields  are  again  approximately  10  V/m; 

2)  The  meter  current,  I , changes  a factor  of  4 with  time; 

3)  The  final  equilibrium  meter  current  is  positive,  not  negative, 
as  seen  in  the  previous  cases. 

(2)  Z1»Z2 

Let  the  metal  electrode  be  beryllium  or  carbon,  and  let  the  dielectric  mate- 
rial have  a relatively  high  effective  Z such  as  Saran  or  polyvinyl- chloride  or 
SiOg  or  AlgOj.  Dielectrics  with  higher  effective  Z can  be  found  such  as  tanta- 
lum oxide,  but  we  have  not  considered  these. 

(a)  200  keV  photons  incident  on  PVC  dielectric  with  beryllium  center 
14  2 

electrode.  101  photons/rn  -sec  are  incident  at  x=  0.  The  curve  fits  to  the 
available  data  are  not  to  be  taken  as  true  curve  fits..  The  available  data  is  not 
complete  but  is  "good  enough"  to  give  us  approximately  correct  fits.  Therefore, 

12.  Frederickson,  A.  R.  and  Burke,  E.  A.  (1971)  IEEE  Transactions  on  Nuclear 

Science  NS18: 162-9. 

13.  Long,  David  M.,  Chadsey,  William  L.,  and  Benedict,  Robert  V.  (1971)  Predic- 

tion of  Dose  Gradients  and  Their  Effects  on  Semiconductors,  AFCRL-V 1-0584. 
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the  results  for  J(x,  t),  E(x,t),  and  I(t)  are  only  indicative  of  the  effects  to  be 
seen.  Exact  predictions  for  200  keV  photon  beams  cannot  be  made  at  this  time 
due  to  the  lack  of  accurate  inte:  * ice  photo- compton  current  and  interface  dose 
data  with  this  material  combination.  Results  are  shown  in  Figure  10. 

PVC  BERYLLIUM  PVC 


X DEPTH  IN  METERS 


BERYLLIUM 


X OEPTH  IN  METERS 


BERYLLIUM 


DEPTH  (METERS) 


TIME  I (SECONDS) 

Figure  10,  0.2  MeV  Photons,  10^  Photons/m^ 
-sec  on  PVC  With  Beryllium  Center  Electrode 


In  this  case,  I(t)  changes  about  8 percent  and  the  electric  fields  peak  at 
—10  V/m.  These  effects  are  significantly  reduced  compared  to  Figure  9.  In 
general,  the  effects  can  be  expected  to  be  a minimum  as  the  materials  approach 
each  other  in  atomic  number.  Strong  effects  could  still  be  seen  for  low  Z - 
low  Z materials  if  the  photon  energy  were  below  20keV.  We  do  not  have  adequate 
dose-depth  data  and  photo- compton  current  data  at  the  lower  energies  to  do  the 
proper  analysis. 

(b)  1.25  MeV  photons  would  produce  effects  qualitatively  the  same  as  the 
200  keV  photons.  We  will  not  include  results  for  the  1.25  MeV  photons  for  this 
geometry. 

7.3  Outer  Electrodes  of  Differing  Atomic  Number 
(1)  Zj  » Zg  as  shown  in  Figure  11. 


Figure  11.  Apparatus  With  Outer  Electrodes  of  Different 
Material 

(a)  1.25  MeV  photons  on,  sav,  Z.  = 29  and  Z„  = 5 (ethylene).  We  can 

• 14  1 . 2 ^ 

obtain  D(x'  and  .!(x)  caused  by  10  photons/m  -sec  at  x=  0.  Again  the  dielec- 

3 

trie  density  is  =*  1.0  g/cm  . The  results,  shown  below  in  Figure  12,  indicate 
that  I(t)  varies  less  than  in  Figure  8.  This  probably  is  due  to  the  mutual  cancel- 
lation of  opposing  effects  in  the  two  dielectrics.  The  maximum  E(x,t)  is  analo- 
gous to  Figure  8 and  occurs  near  the  outer  electrodes. 

14  2 

(b)  200 keV  photons,  10  photons /m  -sec  incident  on  Zj  = 29  and 

Z g =»  5 (ethylene).  The  resulting  .!(x,t),  E(x,t),  D(x),  and  I(t)  are  given  in  Fig- 
ure 13.  The  results  show  again  that  the  time  dependent  fields  and  currents 


30 


produce  a time  dependent  meter  current  I(t),  and  that  peak  electric  fields  of 

O 

nearly  10  V/m  can  be  reached. 


DEPTH  (METERS) 


Figure  12.  1.25  MeV  Photons,  10*“*  Photons/m^ 

-sec  on  Copper  Surrounding  Ethylene 
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TIME  (SKEWS) 


Figure  li.  200keV  Photons,  1014  Pbotons/m^ 
-sec  on  Copper  Surrounding  Ethylene 
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7.4  Further  Cam 

It  is  possible  to  permute  the  material  atomic  numbers  in  several 
additional  ways.  For  example,  we  might  enhance  the  maximum  electric  field  by 
the  following  structure: 


12 

It  has  been  shown  that  this  type  of  material  combination  produces  greater  dose 
gradients  and  a greater  divergence  of  photo-compton  currents  than  the  cases  al- 
ready discussed.  The  increased  dose  gradients  and  current  divergence  will  usu- 
ally result  in  increased  electric  fields  and  significantly  differing  meter  currents, 
I(t).  However,  it  is  not  qualitatively  instructive  to  make  more  permutations,  so 
we  have  not  done  so  here. 

7.S  Mixed  Photon  Spectra 

Interesting  cases  might  occur  when  several  photon  energies  are  simultan- 
eously applied.  The  only  mixed  spectra  case  we  have  considered  to  date  is  simi- 
lar to  Sec.  7. 2(1),  but  with  both  200  keV  and  1250  keV  photons  of  equal  number  in- 
tensity. Many  other  spectra  combinations  might  be  of  interest  but  they  will  not 
be  considered  here. 

In  this  case  we  have  a high  Z center  electrode  surrounded  by  low  Z dielec- 
trics and  low  Z outer  electrodes.  The  total  dose  rate  and  the  total  photo-compton 
current  are  each  the  sum  of  the  200  keV  and  1250keV  dose  rate  and  photo- 
compton  currents  respectively.  For  the  sake  of  variety,  we  have  changed  the 
radiation-induced  conductivity  constant.  Up  to  Up  = 10  sec/ohm-m-rad . (In 
previous  examples.  Up  = 10  ) If  we  compare  the  results  with  Sec.  7.2,  we 

see  that  the  time-to- reach-equilibrium  and  the  electric  field  strength  are  a 
factor  10  lower  in  these  latest  results  due  to  the  change  in  Up  value.  The  re- 
sults are  shown  in  Figure  14. 

Overall,  the  lesults  are  fairly  similar  to  Sec.  7.2;  however,  one  interesting 
result  stands  out.  Notice  that  the  electric  field  changes  sign  twice  as  we  vary 
depth  in  the  right  dielectric.  There  are  two  depths  in  the  dielectric  at  which  the 
electric  field  becomes  zero.  In  all  previous  cases,  there  is  only  one  depth  at 
which  the  electric  field  is  zero.  If  we  look  at  the  plot  of  J(x,  0)  and  consider  its 
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divergence  we  see  that  alternate  layers  of  positive  and  negative  charge  are  being 
deposited  in  the  dielectric  at  various  depths.  If  we  were  to  change  the  right  con- 
ductor to  a high  Z material,  the  photo-compton  current  would  produce  a third 
depth  at  which  the  electric  field  is  zero.  Furthermore,  by  adding  layers  of  dif- 
fering Z dielectric  and  by  irradiating  from  non-normal  directions,  we  might  find 
it  possible  to  increase  the  number  of  electric  field  nodal  points  indefinitely. 

If  the  dielectric  layers  were  relatively  thin,  it  might  be  that  the  low  energy 
components  of  a mixed  spectrum  would  produce  the  predominant  effects.  For  ex- 
ample,  suppose  we  have  a field  effect  transistor  with  a ~ 10  a dielectric.  As- 
sume a fluence  with  a continuous  spectrum  from  1 MeV  to  OeV.  The  electrons 
above  100 keV  would  be  virtually  unperturbed  by  the  interface.  Lower  energy 
electrons  would  be  more  strongly  affected  by  the  dielectric  metal  interface.  The 
abrupt  change  in  electron  scattering  properties  at  the  metal  dielectric  interface 
could  affect  lower  energy  electric  transport  most  strongly.  Also,  the  abrupt 
change  in  photoelectric  cross  section  would  cause  highly  diverging  photoelectron 
currents  for  photons  below  lOOkeV.  The  net  result  in  that  for  thin  dielectrics, 
the  total  charge  deposition  may  depend  most  heavily  on  the  low  energy  components 
of  the  radiation  fluence.  The  effects  of  mixed  spectra  on  thin  dielectrics  presents 
a whole  new  array  of  possibilities. 


7.6  Dependence  on  1!^ 

Earlier  we  noted  that  the  radiation  induced  conductivity  coefficient.  Op  , is 
not  accurately  determined.  It  is  possible  to  obtain  the  functional  dependence  of 
E(x.  t),  I(t),  and  J(x,t)  on  Op  from  the  equations  given  earlier,  but  we  have  not 
done  so.  Instead,  we  have  used  the  data  of  Sec.  7.2  with  several  values  of  0R  to 
determine  the  effect  of  varying  0R  . In  this  way,  our  data  can  be  extrapolated 
for  other  values  of  Op  . It  turns  out  that  if  we  have  calculations  of  Jq(x,  t), 
Eglx.t),  and  IQ(t)  for  a specific  value  of  0R . say  Op  = Opg,  then  we  can  deter- 
mine J(x,  t),  I(t),  and  E<x,t)  for  any  value  of  0R  as  follows: 


E(x,t:0R) 


!^R0e  L J“R_t 

UR  n Uro 


I(t;  Op) 
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‘0  o 


t . 
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Those  relationships  are  valid  only  when  U^  D » Uq  . Notice  that  the  changes  in 
Up  produce  simple  transformations  of  the  time  and  electric  field  intensity  axis. 
There  is  no  change  in  the  spatial  dependence  for  E or  J and  the  initial  and  final 
values  for  I and  J are  unchanged. 

The  physical  explanation  for  the  Up  dependence  is  simple.  The  initial  photo- 
co  'vpton  current,  J(x),  is  fixed  by  the  geometry  independent  of  Up  . The  result- 
ing apatial  charge  deposition  produces  electric  fields  which  increase  with  time 
until  1)IV  J = 0 and  no  more  charge  is  deposited.  The  electric  field  simply  builds 
up  to  the  level  required  to  produce  the  conduction  currents  necessary  to  create  a 
J such  that  DIV  .1=0.  A change  in  Up  does  not  change  the  spatial  dependence  of 
either  the  charge  deposition  or  the  conduction  currents,  but  on.y  changes  the  inten- 
sity of  the  conduction  currents  for  a given  E . Thus,  if  we  decrease  Up,  the  sys- 
tem reacts  by  letting  the  initial  charge  deposition  continue  longer,  until  E is 
finally  increased  sufficiently  to  cause  the  necessary  equilibrium  conduction  cur- 
rents. Thus  the  spatial  dependence  of  J(x,t)  is  unaffected.  Only  the  time  to  reach 
equilibrium,  and  the  final  electric  fields  are  affected  by  changing  Up. 


& APPLICATION  TO  ELECTRONIC  DEVICES 

In  addition  to  the  FET  discussed  above,  coaxial  cables,  dielectric  compton 
diodes,  capacitors,  and  LSI  circuits  are  all  directly  affected  by  charge  injection 

14 

effects.  Time-dependent  currents  to  conductors  have  been  seen  in  coaxial  cables 
and  dielectric  metal  sandwiches.1^  LSI  and  MOS  devices  exhibit  trapped  charge 
in  the  dielectrics  under  irradiation,  and  the  trapped  charge  distribution  is  often 
dependent  on  an  applied  electric  field.  The  mechanism  discussed  ir.  this  paper 
may  play  an  important  role  in  these  devices. 

Our  ability  to  apply  this  mechanism  to  actual  three-dimensional  devices  is 
presently  limited.  It  is  very  difficult  to  calculate  the  photo- compton  currents  at 
a one-dimensional  interface.  Two-dimensional  interfaces  have  not  yet  been  done 
to  our  knowledge,  but  the  potential  exists;  it  is  only  required  that  a Monte  Cax-lo 
program  be  developed  for  a representative  geometry,  such  as  a coaxial  capacitor 
geometry.  For  cases  with  thin  dielectric  layers,  a computer  program  must  be 
especially  dedicated  to  determining  the  currents  and  charge  depositions  in  that 
layer.  However,  the  currents  in  the  dielectric  layer  depend  on  the  electron  scat- 
tering histories  in  material  outside  the  layer.  It  seems  to  us  that  the  thin  dielec- 
tric problem  will  require  computer  programs  especially  dedicated  to  this  one 


14.  van  Lint.  V.  A.  J.  (1970)  IEEE  Trans.  Nuc.  Sci.  NS  17:2 10-16. 

15.  Sullivan,  W.  H.  and  Ewing,  R.  L.  (1971)  IEEF  Trans.  Nuc.  Sci.  NS18:313  (Fig.  5). 
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geometry  so  that  we  maintain  the  spatial  detail  so  necessary  for  calculations 
inside  the  dielectric. 

Let  us  look  at  a thin  layer  problem.  Unpublished  in-house  data  for 
2 

100  Roentgen-cm  / sec  cobalt  60  photon  beam  on  a low  Z - high  Z interface  de- 
posits charge  in  the  first  10®  A of  low  Z material  at  the  rate  of  2X  10  ^C/cm^-R. 
If  we  do  a linear  extrapolation  in  depth,  then  a 10J  A layer  receives  charge  at  the 

rate  of  2 X 10  *®  C'cm^-R.  For  a total  fluence  of  10®  R (cobalt  60),  we  get  a 

-9  2 

charge  deposition  of  2 X 10  C/cm  . Such  a surface  charge  density  on  a layer 

3 . 4 . 

of  SiOg  of  10  A thickness  produces  a voltage  of  ~ 0.1  V.  If  the  layer  were  10  A 

thick,  then  the  voltage  developed  would  be  10  V.  To  calculate  the  charge  deposi- 
3 t 

tion  in  a 10  A layer  requires  doing  better  than  a linear  extrapolation  in  thickness, 
since  the  plot  of  charge  deposition  vs  depth  has  a steep  slope  for  thickness  less 
than  10®  A.  We  are  currently  unable  to  make  a correct  extrapolation  to  10®  A 
thickness  but  we  hope  to  be  able  to  do  so  in  the  future.  The  charge  buildup  in  the 
irradiated  ~10®  A gate  dielectric  of  MOS  devices  may  be  affected  by  the  mechan- 
ism given  in  this  report.  To  this  author's  knowledge,  this  mechanism  has  not 
been  given  consideration. 

The  unstable  response  of  dielectric  compton  diodes  might  be  explained  by 
this  general  model.  However,  the  compton  diode  geometry  generally  in  use  would 
require  a difficult  three-dimensional,  coupled  electron-photon  transport  code. 

The  planar  slab  compton  diode  is  probably  solvable  with  currently  available  data. 
The  vacuum  compton  diode  avoids  the  dielectric  charging  difficulty  but  has  less 
sensitivity.  The  multiple  high  Z - low  Z slab  vacuum  compton  diode  gives 
factor  of  ten  sensitivity  gain  without  the  dielectric  problem  but  appears  to  be 
limited  to  slab  geometry  and  collimated  beams. 

Radiation-induced  conductivity  coefficients  have  been  determined  by  irradiat- 
ing dielectric  slabs  between  biased  metal  plates  and  measuring  the  bias  current. 
We  have  shown  here  that  tremendous  electric  fields  are  built-up  internal  to  the 
dielectric.  Therefore  the  data  taken  with  biased  metal  plates  may  be  incorrect 
since  the  internal  electric  field  strength  was  not  correctly  measured.  It  may  be 
more  accurate  to  use  a dielectric  slab  in  the  geometry  of  Figure  1 with  a known 
J(x)  profile  to  determine  the  induced  conductivity  from  the  measured  time  de- 
pendence of  I(t). 


♦Also  see  U.  S.  Patent  No.  3,  780,  304,  Charge  Accumulation  Gamma  Radiation 
Detector. 


37 


9.  CONCLUSIONS 


A method  has  been  developed  for  determining  the  time-dependent  currents 
and  electric  fields  in  a dielectric- metal  structure  in  one-dimensional  slab  geo- 
metry. Results  are  given  for  200  keV  and  1250keV  photon  irradiation  on  typical 
materials.  Interesting  results  include  time-dependent  electrode  currents  which 

Q 

reverse  sign,  and  electric  fields  which  reach  10  V/m. 

We  have  shown  that  the  dose  depth  dependence  and  electron  current  depth  de- 
pendence strongly  affect  the  currents  to  metal  elements  (sometimes  called  re- 
placement currents).  An  electric  field  builds  up  with  time,  producing  a time- 
dependent  conduction  current.  When  equilibrium  is  reached,  the  high  energy 
electron  current  and  the  conduction  current  sum  to  a net  current  whose  diver- 
gence is  everywhere  zero.  This  net  current  can  be  very  different  from  the  "no 
field"  currents  normally  used  to  predict  "replacement  currents"  in  irradiated 
devices. 

The  relaxation  time  predicted  here  for  organic  dielectrics  is  virtually  identi- 
cal to  the  measured  relaxation  time  reported  for  dielectrics  in  MOS-FET  devices. 
Furthermore,  this  relaxation  time  could  be  the  basis  for  time-dependent  effects 
seen  in  irradiated  capacitors. 

At  this  time  it  appears  prudent  to  experimentally  check  the  results  presented 
here.  If  the  calculations  are  experimentally  confirmed,  then  it  would  be  of  inter- 
est to  obtain  the  information  necessary  to  quantitatively  predict  the  response  of 
some  devices.  This  prediction  might  show,  for  example,  that  capacitor  responses 
will  always  be  time  dependent  and  radiation  history  dependent.  It  might  also  show 
that  MOS  devices  will  always  have  dielectric  charging,  but  th«>t  the  charging  can 
be  minimized  by  a judicious  choice  of  surrounding  material.  To  predict  such  de- 
vice response,  we  require  information  on  7(x)  and  D(x)  in  the  appropriate  geo- 
metry. This  information  is  currently  lacking. 
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Glossary 


V voltage  in  units  of  volts. 

2 

a surface  charge  density  in  units  of  coulombs  /meter  . 

2 

I meter  crrent  in  units  of  amperes/meter  . 

2 

J current  density  in  units  of  amperes  /meter  . 

3 

p volume  charge  density  in  units  of  coulombs /meter  . 

E electric  field  intensity  in  units  of  volts /meter, 

e permittivity  in  units  of  farads /meter. 

A amperes. 

t$0  conductivity  in  units  of  ohm  1 meter  *. 

coefficient  of  radiation  induced  conductivity  in  units  of 
seconds/ohm-meter-rad. 

D radiation  energy  deposition  rate  in  units  of  rads /second. 
(One  rad  = 100  erg/gram). 
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